Temperature evolution of infrared- and Raman-active phonons in graphite 
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We perform a comparative experimental and theoretical study of the temperature dependence 
up to 700 K of the frequency and linewidths of the graphite Eiu and E2g optical phonons (~ 1590 
and 1580 cm~^) by infra-red (IR) and Raman spectroscopy. Despite their similar character, the 
temperature dependence of the two modes is quite different, e.g., the frequency shift of the IR-active 
Eiu mode is almost twice as big as that of the Raman-active E2g mode. Ab initio calculations of the 
anharmonic properties are in remarkable agreement with measurements and explain the observed 
behavior. 
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Thermal properties of sp^ carbon systems such as 
graphene [ll, Q and carbon nanotubes [1, |3] have at- 
tracted significant attention, both systems being excel- 
lent thermal conductors. In these materials, thermal 
transport is dominated by lattice vibrations (phonons) 
and transport can be described properly only once 
phonon scattering mechanisms are taken into account. 
Phonon scattering also plays a major role in non- 
equilibrium phenomena in which nonthermal phonon 
populations of carbon materials are induced by optical 
excitation 0| or during electronic transport 0, la] ■ In 
general, our understanding of thermal properties is heav- 
ily based on theoretical modeling [9l-lll| and case studies 
allowing to validate the models are of primary impor- 
tance. 

The frequency and the linewidth of a phonon, as well 
as their temperature dependence, are measurable quanti- 
ties which provide key information about the interatomic 
interactions [T5|. Indeed, even in a perfect crystal, the 
phonon linewidth is not zero and the phonon frequen- 
cies depend on temperature because of anharmonic terms 
in the interatomic potential (phonon-phonon interaction) 
and, in the same instance, because of non-adiabatic ef- 
fects (phonon-electron interactions) when the gap is suf- 
ficiently small. Raman spectroscopy certainly has a cen- 
tral role in the study of vibrational properties of sp^ car- 
bon materials. However, this technique probes only a 
very small portion of the vibrational degrees of freedom 
and the use of complementary tools is essential to gain 
further insight. Infrared (IR) spectroscopy is a natural 
choice and it can also be used to study planar sp^ sys- 
tems 

The graphene in-plane anti-phase movement of the two 
unit-cell atoms is a Raman-active zero momentum (q=0) 
mode with symmetry i?2g (Raman G peak). The G fre- 
quency depends on the temperature and the G peak Ra- 
man measurement is currently used as a probe of the lo- 
cal temperature in graphene samples P, In graphite, 
the i?2g mode splits in two modes with symmetry E2g 
and Eiu, depending on the relative motion of the two 
graphene planes of the graphite unit cell The E2g 

one (in-phase vibration of the two planes) is Raman ac- 



tive and is associated with the well known G peak at 
^1582 cm~^. The Eiu one (anti-phase vibration of the 
two planes) is IR active and has a slightly higher fre- 
quency ('-^^1590 cm^^). The anomalously large IR inten- 
sity of this mode has been the topic of a recent work 14 1 . 

Given the similarities between the i?2g and Eiu vibra- 
tions, their properties are expected to be similar. On the 
contrary, their linewidths are quite different: At room 
conditions, the full-width at half maximum (FWHM) 
of the Raman mode is ~13 cm~^, while the FWHM 
of the IR one is ~4 cm~^. Based on calculations [llj 
these width differences have been related to the different 
electron-phonon interaction. Moreover, according to Ref. 



lll |. the Raman mode linewidth decreases by increasing 



the temperature (this behavior is quite unusual), while 
the IR linewidth is predicted to increase. Above room 
temperature, measurements of the temperature depen- 
dence of the linewidth and frequency shift of the IR mode 
are needed for a comparison. 

In this Rapid Communication, we present the measure- 
ment of the line- width and the line-shift of the infra-red 
active mode at 1590 cm"-'^ in highly oriented pyrolitic 
graphite in the temperature range 293-523 K. For a com- 
parison, we also show Raman spectroscopy measurements 
of the G peak in a similar temperature range. The results 
are interpreted by means of ab-initio calculations of the 
anharmonic properties. 

Raman spectra were measured using a Horiba Jobin 
Yvon T64000 spectrometer equipped with an Ar laser 
(514.5 nm). The laser was focused on the sample via 
a Mitutoyo 20 x long distance objective. The incoming 
power was selected in order to prevent sample damage 
and ensure a rather high signal to noise ratio. The in- 
strumental resolution obtained by measuring the emis- 
sion spectrum of neon in the i?2g phonon frequency region 
has been found to be 1.6 cm~^ (FWHM of a Gaussian 
profile). Spectra were recorded by spanning the region 
( 1 200- 1 700 cm^ ^ ) . The sample was mounted inside an in- 
house made oven composed of a cartridge heater tightly 
inserted in a cylindrical copper piece with a little squared 
open cavity for the incoming light. The heater was con- 
nected to a power supply equipped with a thermome- 
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FIG. 1: (Color online) Graphite measured Raman spectra and 
IR optical conductivity spectra. Spectra taken at different 
temperatures (T) are vertically shifted for clarity. Each peak 
is fitted with a single Lorentzian (thin dotted lines). 



ter allowing to control the temperature of the cartridge 
heater. The precision was ±0.1C. 

IR spectra were measured at the AILES infrared beam- 
line Hirl of synchrotron SOLEIL (Saint Aubain, France) 
using a Bruker IFS125 spectrometer working under vac- 
uum. The temperature was set by using a in house made 
oven similar to the one used for the Raman experiments. 
The sample was aligned in reflection geometry allowing to 
acquire reflectivity measurements of the a6-plane trans- 
verse excitations at increasing temperatures. The ref- 
erence was measured at room temperature using a gold 
mirror mounted in such a way as to replace the sample 
inside the oven. The entire working set-up used was as 
follows: synchrotron light as the source, KBr as the beam 
splitter and mercury cadmium telluride (MCT) as the 
detector. With this configuration we acquired spectra in 
the frequency region between 700 and 3500 cm~^ with a 
resolution of 0.1 cm"'^. Reflectivity data were treated by 
Kramers-Kronig transformation to extract the real part 
of the optical conductivity (t{uj). The Raman spectra and 
the conductivity IR spectra, corrected by the respective 
background, were fitted with a single Lorentzian profile 
(Fig. [1]) to obtain the width and position of the peaks. 
For a discussion on the asymmetric shape of the IR peak 
see [11]. 

In a defect-free perfect crystal, phonons have a finite 
lifetime because they can decay by anharmonic phonon- 
phonon (ph — ph) scattering or by electron-phonon (e — 
ph) scattering, this last process being possible in graphite 
as the electronic gap is zero. The intrinsic linewidth of 
a phonon can thus be written as a sum of two terms 




7(r) = jPf'-P'^ (T) (T) (see e.g. and references 
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FIG. 2: (Color online) Upper: Graphite phonon frequency as 
a function of temperature. Symbols are the present measure- 
ments from Raman (triangles) and IR (dots). Solid lines are 
calculations. Lower: Linewidths (full width at half maximum) 
corresponding to the upper panel measurements. Dashed 
lines are the contribution to the linewidth due to anharmonic 
phonon-phonon scattering. Continuous lines are the total in- 
trinsic linewidth (phonon-phonon -I- electron-phonon, see the 
text). The calculated Raman shift and the calculated Raman 
and IR linewidths are taken from Ref. [Tm ■ 



therein), where T is the temperature. The linewidth in a 
real sample can be larger because of the presence of other 
scattering mechanisms. Concerning the temperature de- 
pendence of the phonon pulsation ujph, at the lowest order 
one can distinguish the following contributions [l^ : 

ujph^0J0 + Auj'''{T) + Auj''P{T) + ALj^P{T) + O{h^), (1) 

where ujq is the harmonic pulsation at the equilibrium 
lattice parameters, Auj'-'^{T) accounts for the variation 
of the harmonic pulsation by varying the lattice param- 
eters (as a consequence of the thermal expansion), and 
Aui^P{T) and Auj'^p{T) can be interpreted as due to an- 
harmonic scattering involving, respectively, three or four 
phonons (see e.g. [l^ [13] and references therein). The 
three shifts are proportional to h and are, thus, expected 
to be of the same order of magnitude. Other terms are 

In the present study, we compare the measured 
linewidths with those calculated in [ll[. Reference [ll| 
reports the calculated anharmonic lineshift for the Ra- 
man mode only. Here, we determine the anharmonic 
lineshift of the IR active mode using the same approach 
as Calculations are performed within density 

functional theory (plane-waves and pseudopotential ap- 
proach) with the QUANTUM ESPRESSO package [l|. 
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Phonon dispersions are obtained within the hnear re- 
sponse approach of [Tgj . Anharmonic phonon-phonon 
scattering coefficients are determined thanks to the 2n-|-l 
theorem as implemented in \2d\. All the computational 
details have been described in 11], and, in particular, 
we use the local density approximation and equilibrium 
lattice spacing parameters (a — 2.43A, c/a — 2.725). 
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FIG. 3: (Color online) Calculated lineshift of the E2g (Ra- 
man active) and Eiu (IR active) T phonons of graphite. The 
total lineshift (continuous line) is decomposed into the three- 
phonon, four-phonon and lattice expansion contributions. All 
the contributions are vertically shifted so as to be zero for 
r=OK. The r=OK values are given in the text. Calculations 
for the Raman mode are taken from Ref. [flj . 

Figure[3]compares the measured and calculated param- 
eters for the Raman (i?2g) and IR active (Eiu) modes. 
Measured linewidths are in good agreement with the cal- 
culations of [ll|: The IR linewidth increases when in- 
creasing the temperature while the Raman decreases. 
These trends are directly visible on the raw data of 
Fig. [1] This behavior was interpreted in [ll| consider- 
ing the following: (i) The phonon-phonon contribution to 
the linewidth, 'jP^~p'^^ is very similar for the two modes 
and increases with the temperature (Fig.[T|); and (ii) the 
ryC-ph contribution is significantly more important for the 
E2g mode and provides the overall temperature decrease 
of the linewidth. 



where is a sum on the Brillouin zone, f'^'^" is the unit 
cell energy, n^j is the Bose-statistics occupation of the 
phonon (qj), and d is the Dirac distribution. Equation 
([2]) depends on the temperature only thr oug h the terms 
Hqj. The associated four-phonon shift jl2l [iTj can be 



Figure [2] also shows the line-shift. The theoretical 
curves are vertically shifted to best fit the experimen- 
tal data, obtaining, for T=OK, 1592.6 and 1582.9 cm^^ 
for Eiu and i?2g, respectively. These can be considered 
as the extrapolated T—OK experimental frequencies and 
their difference is 9.7 cm^^. At the equilibrium lattice 
spacing, the (i?2g) frequency is 1614.2 (1604.9) cm~^. 
To compare with the r=OK experimental value, one has 
to add the calculated anharmonic shifts at T=QK, which 
are Aw''=(0) = -17.7 (-17.7) cm"!, Acj^p = -14.2 (-14.1) 
cm-i and Auj^p = -1.3 (+1.3) cm^^ for the Eiu {E2g) 
mode. The resulting frequency difference is 6.6 cm~^, 
which underestimates by 3 cm~^ the experimental 9.7 
cm~^ value. By considering frequencies at the general- 
ized gradient approximation (GGA) level, [2l| the agree- 
ment slightly worsens. 

By increasing the temperature from 293 to 523 K, 
the measured IR frequency shifts by —7.3 cm^^, which 
is more than twice the corresponding Raman shift, 
—3.3 cm^^, in agreement with calculations (Fig. [2]). To 
understand this behavior, we decompose the lineshift of 
the two modes into the three T dependent components 
defined in Eq. [T] (see Fig. The different behavior of 
the two E2g and Eiu modes is almost entirely determined 
by the corresponding 4-phonon scattering terms (Fig. [3]). 

Why are the four-phonon shifts of the two modes so 
different? To answer this question, we need to intro- 
duce some concepts. Let us consider the interatomic 
potential energy 5*°*({usq(R/)}), where Usct(Rz) is the 
displacement from the equilibrium position of the s- 
th atom in the crystal cell identified by the lattice 
vector R/ along the a Cartesian coordinate. For a 
phonon with wavevector q, branch index i, and energy 
fiijjqi, we define the adimensional phonon displacement 

as Uqi = v^2A'/sWqi/?lUsQ(qi)"sa(R0e"*'''^', 

Vsai^i) being the orthogonal phonon eigenmodes normal- 
ized on the unit cell, Mg the atomic mass, and N the 
number of q-points describing the system (or unit cells) . 

For a phonon q=0, and branch index i, we define the 
four-phonon spectral function as: 



(2) 

I 

written as 

^^of = / doj S'of (w). 
Jo 

The fourth order derivatives in Equation ^ can be in- 
terpreted as scattering coefficients among two phonons 
Oi and the two phonons -qj,qj. The spectral function 



I 

'5'of (^) = TTTTr Y] 7^ J. 7^ (2nqj- + 1)5{UJ - Wqj-), 
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FIG. 4: (Color online) (a): Scheme of various F (q=0) phonons of graphite, (b): Graphite phonon dispersion and related 
vibrational density of states (VDOS). (c): Spectral decomposition of the four-phonon component of the frequency shift [as 
defined in Eq. ([2])] of the £29 (Raman) and Eiu (IR) phonons at temperatures and 300 K. 



separates the contributions from phonons virith different 
energies. 

The spectral functions of the two E2g and Eiu modes 
are very similar [Fig. SJc)], the major differences being 
near 115 cm~^, where the Eiu spectral function shows a 
peak absent in the E2g one. This peak strongly depends 
on the temperature (in the temperature range consid- 
ered) and accounts for the difference in the two E2g and 
Eiu shifts. This peak is associated with the maximum of 
the graphite phonon density of states (VDOS in Fig. 01) at 
115 cm"""^, and this maximum is due to the optical out of 
plane ZO' phonon branch. We remind that, in graphene, 
the vibration of the atoms perpendicular to the plane, is 
associated with the ZA acoustic branch. In graphite, this 
branch, splits into a branch in which the two graphene 
planes vibrate in-phase, ZA, and a branch in which the 
two planes vibrate anti-phase, ZO', [Fig. |31Ja)-(b)]. The 
ZA graphite branch is acoustic, while the ZO' is optical 
with a 114 cm^^ frequency at q=0. The four-phonon 
anharmonic coupling [the fourth order derivative in Eq. 
([2])] between the mode and the ZO' branch is much 
stronger than that between the E2g and the ZO'. The 
presence (absence) of this coupling explains the presence 
(absence) of the 115 cm~^ peak in the spectral decom- 
position of Fig. I3J; for the £'i„(i?2g) mode. 

We remark that graphene bilayer also presents an IR 
active Eiu optical mode and a ZO' mode with a fre- 
quency of ~ 80 cm~^ [22I, smaller than in graphite. A 
smaller frequency is associated with a stronger temper- 



ature dependence of the ZO' phonon occupation (with 
respect to the graphite ZO'). This fact, suggests that 
the IR active Eiu mode of the bilayer should present 
a temperature dependence of the shift stronger than in 
graphite. 

To conclude, we measured the linewidth and frequency 
shift of the F optical phonons Eiu (IR- active) and E2g 
(Raman active) of graphite as a function of temperature. 
Despite the similarities between the two Eiu and E2g 
vibrations, their temperature dependence is quite differ- 
ent. The Eiu linewidth is almost three times smaller 
than the E2g one (at room temperature) and, while the 
Eiu linewidth increases by increasing the temperature, 
the E2g one slightly decreases. Both behaviors are in 
good agreement with the ab initio calculations of [Tl| . 
The shifts of the two modes are also very different: By 
increasing the temperature from 293 to 523 K, the Eiu 
mode shifts by more than twice the corresponding E2g 
value, in agreement with the present calculations. This 
difference is explained by the presence (absence) of the 
four-phonon anharmonic coupling between the Eiu (£-29) 
phonon and the ZO' phonon branch at ~ 115 cm~^. 
These findings confirm the accuracy and predictive power 
of density functional theory calculations in determining 
anharmonic phonon-phonon scattering in sp^ carbon sys- 
tems, a key step in developing realistic models for thermal 
transport Q . 

Part of the calculations were done at IDRIS, project 
096128. 
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